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Abstract 

We employ the Weizsacker- Williams method of virtual quanta to 
study disintegration of nuclei in magnetic field. We explore a variety of 
field configurations and conclude that for the energy range of interest 
for applications to cosmic rays (10 18 — 10 21 eV) such disintegrations 
are not a significant source of energy or flux loss for any realistic 
acceleration mechanism. 

1 Introduction 

Understanding the interactions of cosmic rays with cosmic backgrounds is 
essential for solving the puzzle of ultrahigh-energy cosmic rays 0, as well 
as for understanding the most powerful natural accelerators, such as active 
galactic nuclei, gamma-ray bursts, etc. Since their discovery, cosmic rays with 
energies beyond the Greisen-Zatsepin-Kuzmin cutoff 0, E > 10 19 6 eV, have 
presented astrophysics with several important and unresolved problems. The 
composition of ultrahigh-energy cosmic rays remains unknown. In particular, 
they may be protons, photons, or nuclei. In the latter case, interactions with 
the background photons cause nuclear disintagration B. 

A number of possible acceleration mechanisms have been proposed. Some 
of them, such as neutron stars [|J], have very large magnetic fields. To a 
highly relativistic particle, such magnetic fields are much larger and they are 
accompanied by a matching electric field. Thus, the magnetic field acts like 
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an intense photon field. It is therefore natural to ask whether such magnetic 
fields may cause energy losses in relativistic nuclei, through photodisintegra- 
tion. 

In this paper we explore the possibility that the static magnetic fields 
of an accelerator may cause photodisintegration of relativistic nuclei. We 
will use the Weizsacker- Williams method of virtual quanta, as described in 
reference ||. The same method has been used to study photodisintegration 
of nuclei in the Coulomb elctric fields of target nuclei 0. In sections 2-4 we 
set up the method by first transforming the magnetic field to the rest frame 
of the accelerated particle, second, decomposing this field as a collection of 
virtual photons, and third, studying the effect of this spectrum of virtual 
photons on the propagating nuclei. In section 5 we use our method to find 
the effective photon spectrum of a variety of magnetic fields. In section 6 we 
use some simple cross section data to place limits on the required magnetic 
fields to get significant photodisintegration. In section 7 we generalize the 
method to apply to accelerating particles. Finally in section 8 we conclude. 

2 Transformation of the Magnetic Field 

If we attempt to study the disintegration of a relativistic nucleus in a static 
magnetic field from the rest frame of the field, we run into a number of 
problems. In this frame the field does not look like a collection of quasi-real 
photons and therefore the experimental data on photodisintegration would 
be useless to us. We would then have to resort to a full field theoretic 
calculation of the nuclear disintegration. Such theory is both complicated 
and incomplete. Therefore, this path is not very helpful. 

Instead we will begin by transforming our magnetic field to the rest frame 
of the nucleus. In this frame, the field looks very much like a collection of 
real photons (hence the term quasi- real). We express the field as a collection 
of quasi-real photons and use the known experimental data on photodisin- 
tegration to calculate the likelihood that the nucleus will disintegrate in the 
magnetic field. 

We begin with a static magnetic field B(x, y, z). For simplicity we will 
Lorentz boost this in the x direction, so that the argument becomes x — > 
7(x + vt), y — > y, z — > z. The transformations of the fields are also simple 
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and can be found in any standard text such as M. 



v 



E' ± = 7 E1 + -XBJ (1) 



c 

Ej| = E|| (2) 

B' ± = t (b_ l -^xE_ l ) (3) 

Bj, = B,| (4) 

We now substitute this expression into our original form for the magnetic 
field and find 

B' = 7 B± + B|| (5) 

E' = 7- x B. (6) 

c 

Here we note two important aspects of this electromagnetic field. First, 
the parallel component of the magnetic field is irrelevant. It has no factor of 
7 enhancement, and it has no partner electric field and therefore has no effect 
on a particle at rest. This is expected as magnetic fields along the direction 
of a charged particle's motion generally have little or no effect. 

Next, and more importantly, we notice that the electric field and the 
perpendicular component of the magnetic field are perpendicular to each 
other and, for ultrarelativistic particles (v ~ c), approximately the same 
magnitude. Therefore, they can be viewed as a superposition of quasi-real 
photons. This is the subject of our next section. 



3 The Power Spectrum 

We now proceed to calculating the power and number spectrum of the quasi- 
real photons. Using eq. (15.51) of Ref. ||, we find that 

Where is the Fourier transform of the electric field defined by 

E(w) = -=/ E(t)e iUt dt. (8) 



'2-k 
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Substituting in equation |B| for E(£) and changing variables we find 

/oo 
p x B( 7 (x + Pct),y, z)e iwt dt 
-oo 

1 -U)X f°° -wx 1 

—=e^ P x B(x',y, z)e~ c dx . (9) 

\/27T J-oo 



PcV2n 
Where = * 

This equation was derived with the assumption that the particle follows 
a straight path through the magnetic field, but it is easy to generalize it. 
The integral along x will simply become a path integral along the classical 
trajectory of the particle r(£), where I is the path length. Then, droping 
factors of P and the irrelevant phase factor in front of the integral, we get 

E(w) = — ^= I px B(r (i^e^di. (10) 

CV 27T ./path 

4 Photon flux and disintegration probability 

It is now straightforward to calculate the disintegration probability for the 
nucleus. First, we relate the intensity to the flux of photons using 

dl dn 
du du 

Theerefore, 

dn 1 c . , 9 
— = Ew 

duj hLu2n l v n 

Next we multiply this flux by the experimentally determined cross section for 
photodisintegration of the nucleus to find the differential decay probability. 

1 c 

d-Pdecay = a^(u)dn(u) = cr y (u) — — I E(^) 1 2 duj. (11) 

llUJ Z7T 

One change should be made at this point. So far we have neglected 
the possibility that the particle has already decayed. Thus, if we integrate 
equation ([TT]) we may get answers greater than 1, which is impossible. To 
rectify this we need to multiply the right hand side of the equation by the 
probability that the particle has not disintegrated. Thus, equation (Till ) would 
be better written as 

dPdccay dP SUTV { ve -P surv i vc (T^, (cj)rf'Ti(co') . 
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This can then be integrated to find that in fact the correct relation is 

In P surviv0 = - / a^(u)dn(uj) = - aJu) — — |E(c^)| 2 da;. (12) 
Jo Jo nu 2tc 

5 Various field configurations 

We now examine a variety of field configurations to see how E(u;) depends 
on B(t). We define k = ^. 

5.1 Gausian B-field 

Let us suppose that the field is of the form B = j£> exp j— ^-j and the nu- 
cleus is traveling at ultrareletivistic speed in the x direction. From equation 
flllj) we find that 

E{u) = tJfa^{-M} (13) 

We note that for sufficiently small values of a this is relatively fiat, but that, 
for large values of a, E(u) falls off very quickly with increasing uj. However, 
for large enough values of gamma even very large values of to may have 
noticable contributions. 



5.2 Constant 5-field with various cutoffs 

Next we consider a -B-field that is relatively constant over a large region of 
size a and which falls to over a much smaller region of size b. The method 
we will use is to convolve a constant magnetic field with some approximation 
to a delta funtion, thus introducing a smoother cutoff. This provides a nice 
way to compare different kinds of cutoffs on the same distance scale. Also 
the Fourier transforms are then fairly easy to perform because the Fourier 
transform of a convolution is just the product of the Fourier transforms. For 
the rest of this subsection we will leave off the factor of l/c^/2n in equation 

5.2.1 Constant Magnetic Field 

We'll start with a constant magnetic field: 

' B x G (-a/2, a/2) 



B (a;) 



otherwise. 
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The Fourier transform of this is simply 

B / ka\ 
B(*) = f 2«n^ T J. 

If ka <C 2ix then this is nearly constant, however I expect that is unlikely. In 
the opposite limit ka 2tt the oscilations in the sine function are enough 
that we can just replace it by it's root mean square value of \/2. 

It is worth noting that for this configuration the -B-field itself is discon- 
tinuous and the field falls off as 1/k for large values of k. 



5.2.2 Square pulse 5-function 

From now on we will modify B by convolving it with some function g(x, 6) 
that has a characteristic width of b and becomes a 5-function in the limit 
that 6^0. (Basically this is just requiring that the area be normalized to 
1.) If g were a true 5- function this would just give the same form we found 
before. To state our assumptions mathematically 



B(x) 



/oo 
g(x - x')B (x')dx' 
-oo 



J — I 



g{x - x') = 1. 

As our first approximate 5-function we will use a square pulse. 



g(x) = 



\ xe (-6/2, 6/2) 
otherwise. 



The Fourier transform is simply g(k) = SUi ^l^ ) so that the Fourier transform 
of the B-field is 

B ka sin (kb/2) , , 

B w=f 2sm Yi^r- < 14 » 

Note that the small k ik <^ 1/6) behavior is essentially unchanged, but 
that for large values of k (k 3> 1/6) the field now falls off like 1/k 2 . Also 
note that whereas before the B-field itself was discontinuous, now only its 
first derivative is discontinuous. 
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5.2.3 Triangle pulse 



The next step is to consider a continuous approximate 5-function with a 
discontinuous first derivative. This will make the B-field continuous through 
the first derivative. We find the following results. 



g{x) = < 



1(1 + 1) 

Kl "I) 




x e (-6,0) 
x e (0,6) 
otherwise. 



2(l-cosA;6) fsin(kb/2)^ 




(kb/2)\ 



(»/2) ) (15> 
Again the behavior for small kb is relatively stable, while for large kb the 
function falls off as 1/k 3 



5.2.4 Parabolic pulse 

Now we make the approximate (^-function from parabolas joined together in 
such a way that the function and first derivative are both continuous. Thus 
the B-field will have continuous 1st and 2nd derivatives. Thus, 

UW-* 2 ) a; G (-6/2, 6/2) 

M X ~W (6/2,36/2) 

_i_( x + | 6 )2 xe (-36/2,-6/2) 

otherwise. 

Clearly, the observed pattern is continuing. I have chosen a very par- 
ticular form for the parabolic pulse. It is the pulse derived by convoluting 
the square and triangular pulses, and thus gives the very simple form listed. 
However, other parabolic pulses of width 6 and unit area share the property 
that the Fourier transform goes smoothly to 1 for small k and falls off as 



g{x) = < 
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fc~ 4 for large k. I am uncertain how to prove this pattern in general, and the 
specifics are getting increasingly more difficult so I will not attempt smoother 
cutoffs of this form. 



5.2.5 Gaussian Pulse 

It is; however, simple enough to explore the behavior of a particular perfectly 
smooth cutoff by using a gausian pulse as our approximate 5-function. This 
will yield a magnetic field that falls off quickly, but with no discontinuous 
derivatives. We expect to find that B(/c) falls off more quickly than any 
power of k. This is exactly what we find. 

g{k) = exp — 

B . ka k 2 b 2 
B {k) = -£-2 sin — exp — (17) 



5.2.6 Summary and Discussion 

Table [I] summarises B(/c) for the various field configurations discussed in this 
subsection. The last collumn simplifies the form of B(&) by substituting a 
root mean square value for the oscillatory part. 

We should be clear at this point about what is meant by a discontinu- 
ity. Certainly in any real situation we would expect the field and all of its 
derivatives to be smooth down to the atomic scale. However, from the point 
of view of the Fourier transform, any change that occurs over a scale that is 
small compared to the wavelength, A = 2ir/k, will have the same effect as 
a discontinuity. This can be seen in table |I[ Each of the forms reduces to 
the case where the field itself is discontinuous if the smoothing takes place 
over a distance that is small compared to the wavelength. In other words if 
kb <C 1 or equivalently b <C \/2ir then the smoothing does not matter and 
B(ife) = &*2sin*f. 

By virtue of the Fourier transform, the more severe the discontinuity in 
B(/c), the slower the fall-off in B(/c). If the first discontinuity in B or it's 
derivatives appears at order or the zth derivative of B, then for large k 
the B(/c) appears to fall off as k~( t+1 \ 
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Table 1: B (k) 



Smoothing Pulse 


First Discontinuity 


B(k) 


simplified B(k) 


True 5-function 


B 


f> 2 sin !f 


>/2B„ 

k 


Square Pulse 


B' 


B Q„; n feasin(fc6/2) 
k 2 (kb/2) 


2B 

k 2 b 


Triangle Pulse 


B" 


Boo ka ( sin {kb/2) \ 2 
k Z blIi 2 V (kb/2) J 


k 3 b 2 


Parabola Pulse 


B'" 


Boonjv, ka /sin(fc6/2)\ 3 
k Z blIi 2 V ikb/2) ) 


2 v / 10Bo 
fc 4 fe 3 


Gaussian Pulse 


None 


^2 sin k ? a exp 


V2B n 

k C ^P 4 



The pattern that we see here appears to be a general property of the 
Fourier transform that we are taking. The Fourier transform is most sensitive 
to structure (changes in magnetic field, discontinuities in the field or its 
derivatives, etc.) that are on the same scale or smaller than the wavelength. 
The smoother the field is on this length scale, the more heavily surpressed the 
Fourier transform is. This means that in order to see a significant effect on 
nuclei, the magnetic field will need to have significant structure on the scale 
l^dis where is a typical wavelength for photo- disintegration. As we shall 
see later, for the highest known cosmic-ray energies, this requires structure 
on the order of a few milimeters. However, the length scale increases with 
increasing cosmic ray energy. Also, repeating variations in magnetic field, 
which would result in a smaller average field for other measurements, would 
increase the size of B(/c) 

5.3 Dipole Field 

A dipole field can easily be represented by a scalar potential. This is a bit 
unconventional for a magnetic field, but it can be done, so long as there are 
no currents. (This way V x B = 0.) The field B is then simply B = V</>*. 
We can now do the Fourier transform on the single potential rather than on 
the three components of the field. To justify this, one can easily show that if 

B(x,y,z) = V<f>*(x,y,z), 

then 

B Xjy (x,y,k) = d XiV 4>*(x,y,k) 
B z (x,y,h) = ik(j)*(x,y,k). 



9 



The last term will not affect our calculation because it has no corresponding 
E-field. 

For a dipole field exterior to all currents generating the field, the potential 

is 



(18) 



or in terms of components 



p x x + p y y + p z z 
[z 2 + p 2 ]i 



Where I have used p = \/x 2 + y 2 for simplicity when appropriate. The 
Fourier transform of this configuration is 



<P*(x,y,k) = 2 



k 

(p x x + p y y)-K 1 (kp) - ip z kK (kp) 



(19) 



Where Ki is a modified bessel function. 

For kp 1, the following assymptotic expansion is valid: 



KAz) 



' At 2 - 1 
1 + — ; + ••• 



Therefore, one can simplify eq. (p~9|) : 

fT± ■ p 



(j)*(x,y,k) 



P 



ip 2 



<2ixk 



-kp 



(20) 



We can now take the gradiant of this (in the x, y-plane) to find B 

B( P ,e,k)-- 



i2itk 



-kp 



P 



1 



kiipz — yU± cos 8) H (ip z — 3p± cost 

2p 



+ 



p± 



It is worth noting that the leading behavior in p for k ^> 1/p is 
B(A;) = pk{ip z — p± cosi 



j 2nk 



-kp 



P 



(21) 



(22) 
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6 Required B(fc) 



We now estimate the required 5-field to get significant photodisintegration 
due to the static fields, for the various field configurations that we have de- 
scribed. We use the approximate cross sections from Ref. |J. The cross sec- 
tions are given in units of the classical dipole sum rule or J2d = 59.8^ MeV- 
mb, and for almost all of the species listed, the integrated cross section is 
about 1 J2d f° r single nucleon emission and an additional (.1 — .2) J2d f° r dou- 
ble nucleon emission. The cross sections generally reach a maximum around 
20 MeV with a width of 8 — 12 MeV. Clearly, the behavior is similar for all 
of these nuclei. Therefore, we will simply study an extremely simplified case. 
We will take J2d = 60^p and use N, Z, and A for iron. We will take the in- 
tegrated cross section to be 1 J2d centered at 20 MeV with a width of 8 MeV. 
We will ignore any variation in B(fc) over this range and simply substitute 
the value at 20MeV. All of these approximations should be correct to within 
20% for iron and will not result in any qualitative changes in the result. 

Using this data I calculated the magnetic field necessary to get a decay 
probability of 1 — e _1 for an iron nucleus with a gamma factor of 10 10 or 
an energy of about 5 x 10 20 eV. The results can be found in table |2|. These 
results can be modified for other values of gamma by multiplying k times 
We find that k is approximately 100/cm. (Incidentally, this is close to 
the shorter wavelength part of the CMB, which also starts to be important at 
these energies. ||) Entering this in equation ( |I2"D and restoring the l/(cy/2n) 
we find that the necessary B{k) is 

B(Jfe)Jfe = 1.7 x 10 7 gauss. 

Substituting this into our results from table [l] we can make estimates for 
the required magnitude of the magnetic field for each of the field configura- 
tions listed. (Note: Discontinuity refers to the first derivative that changes 
abruptly on the scale of 1/k.) 

It is difficult to imagine objects which would contain fields of this magni- 
tude. Since one would expect the magnetic field strengths to be continuous 
and to fall off over scales of at least a few km, we need magnetic fields of at 
least a 10 14 gauss. For smoother decreases, larger objects, or motion that is 
not perpendicular to the magnetic field the prospects get worse. However, 
for higher energies or magnetic fields with repeated large variations the effect 
may become noticable. 
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Table 2: Required B-field for photodisintegration 



Discontinuity 


Required B field 


Value 


B 


B(k)k 


1.2 x 10 7 gauss 


B' 


B ^ k kb 


8.5 x 10 13 ^6 

km 


B" 




4.9 x io 2 °sp^6 2 


B'" 




2.6 x 10 27 ^6 3 

km 


None 


^ exp( ( T) 


1.2 x 10 7 exp( 2 - 5x k 10 2 13fc2 ) gauss 



7 Accelerating Particles 

All of the above treatment assumes that the particles are moving at a 
constant speed through the magnetic field, and that photodisintegration is 
caused by changes in the magnetic field. It is also possible that a particle is 
accelerated in the presence of a magnetic field, and that the changes in its 
perceived magnetic field due to this acceleration could cause photodisinte- 
gration. We discuss this possibility here. 

If a particle feels a continuous proper acceleration a, (or equivalently, if 
it is pushed with a constant exterior force of F — mod) then the relation 
between 7, x, and the proper time r is as follows. 

7(t) = cosh (ar + b) 

x(t) = - cosh (ar + 6) + C 
a 

Here b and C are integration constants determined by the innitial conditions, 
and we have set the speed of light equal to I. 

Hyperbolic functions can be difficult to work with, however, as long 
as cosh (ar + 6) ^> 1 we can make the substitution cosh (ar + 6) ~ 
I exp (ar + 6). Then the above equations simplify to 

i(t) = 7 exp (ar) (23) 

x(t) = —[exp (ar) - 1] + x . (24) 
a 

We would like to follow the same procedure as in sections 2 and 3 to 
find a power spectrum for the effective electromagnetic field that the particle 
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sees. Most of this except for the Fourier transform will proceed as before. 
However in finding 



BM 



BW 



a number of changes are in order. First, during acceleration the relation- 
ship between x and r is now described by equation (121). Second, while it is 
reasonable to approximate a particle as always moving at a constant veloc- 
ity, it is not reasonable to approximate it as always accelerating. Therefore, 
we must consider the particle's motion before, during, and after the accel- 
eration. It may be tempting to simply assume (correctly) that most of the 
contribution to the Fourier transform will come from the acceleration, and 
therefore ignore the incoming and outgoing motion. However, to do so would 
introduce spurious contributions as we shall see. We will therefore assume 
that the particle comes in from — oo to point xq with 70, that it accelerates 
between xq and X\ from 70 to 71, and that it finally leaves to 00 with 71. The 
correct form for the Fourier transform, after changing integration variables 
is then 



(25) 
(26) 
(27) 



B(u;) = / B{x)e i ^ {x - xo) dx 

J — CO 

+ / B(ar)[— (x-x)Y*dx 



+ 



'11: 

'To' 



i — lx—xi) 



dx. 



_ 21 

a 



Where x = x — = x\ 

If we assume that the magnetic field is constant except that it falls off 
slowly at ±00 so that we can ignore the integration limits there, we find that 



BM 




1 



iu + a iu 




oj* — iau 



■ (28) 



Where the first term in [ ] comes from the acceleration integral and the sec- 
ond term comes from the incoming and outgoing integrals. Notice that as 
mentioned before, the overall behavior in uj is softer than in either term 
alone. If we had integrated only over the acceleration region we would have 
essentially introduced a fictitious abrupt cutoff in the magnetic field. The 
incoming and outgoing motion cancels this effect. We can see that if the 
acceleration is much larger than uj the falloff with uj is fairly slow, and the 
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problem becomes similar to an abrupt falloff in the B-field. As can be seen 
in table [| this may be an achievable magnetic field in some extreme environ- 
ments. However, such high accelerations require that the particle's energy 
change by nine orders of magnitude in a distance of a few microns. Thus, 
we expect much softer behavior, which would require much larger magnetic 
fields to achieve photodisintegration. 



8 Conclusion 

It apears unlikely that static magnetic fields would cause significant photo- 
disintegration of nuclei, unless the fields are rapidly changing or the particles 
are accelerating, on scales of a few microns to a few millimeters. Not coin- 
cidentally, this is the same scale as the shorter wavelengths of the CMB || 
which cause attenuation of such high energy cosmic rays as they traverse the 
universe. 

While this suggests that static magnetic fields are unlike to cause large 
energy losses in cosmic ray nuclei, the effect may become more significant at 
higher energies, or if the magnetic fields experience repeated variations. 
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